Numerical study of the local contractivity 
of the <I>Q mapping 

by 

Marietta Manolessou 
and 
Sofiane Tafat 

E.I.S.T.I. 

Avenue du Pare 95011 Cergy-Pontoise-Cedex 
December 18, 2012 



Abstract 

Previous results on the non trivial solution of the $^-equations of 
motion for the Green's functions in the Euclidean space (of < r < 4 
dimensions) in the Wightman Quantum Field theory framework, are 
reviewed in the 0— dimensional case from the following two aspects: 

• (cf-|l]) The structure of the subset $ C S characterized by the 
bounds signs and "splitting" (factorization properties) is refhned 
and more explictly described in terms of a new closed subset 
$0 C $ C S. Using a new norm we establish the local contrac- 
tivity of the corresponding <i>Q mapping in the neighborhood of a 
nontrivial sequence Hq G $o- 

• A new $g iteration is defined in the neighborhood of the sequence 
Ho e $0- 

In this paper we present the results of our numerical study, so: 

— a) the stability of $o *• e- the splitting, bounds and sign prop- 
erties is clearly illustrated in the neighborhood of Hq € $o • 

— b) the rapid convergence of this iteration to the fixed point 
is perfectly realized thanks to the new starting points of the 
iteration. 
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1 Introduction 



1 A new non perturbative method - the new recent results 

Several years ago we started a program for the construction of a non trivial 
$4 model consistent with the general principles of a Wightman Quantum 
Field Theory (Q.F.T.) [Ij. In reference [2] we have introduced a non pertur- 
bative method for the construction of a non trivial solution of the system of 
the equations of motion for the Green's functions, in the Euclidean space 
of zero, one and two dimensions. In reference [3] we tried to apply an exten- 
sion of this method to the case of four-(and a fortiori of three-)dimensional 
Euclidean momentum space. 

The general aspects of the method together with its comparison and 
validity arguments with respect to other non perturbative methods are pre- 
sented in these previous references, and in particular in the theoretical aspect 
of this last study [4J. 

In this paper we present the numerical study of a new <I>Q-Iteration. 

We use the zero dimensional analog of the system of equations of mo- 
tion introduced in the previous papers. The reasons that motivated us for 
a study in smaller dimensions and not directly in four, were the absence of 
the difficulties due to the renormalization and the pure combinatorial char- 
acter of the problem in zero dimensions. Another useful aspect of the zero 
dimensional case is the fact that it provides a direct way to test numerically 
the validity of the method. 

What are the new developments in the present zero dimensional study: 

1. (cf.[4|)The "new $o subset is explicitly given in terms of the "splitting" 
sequences upper {5n,maxi^)} and lower {5n,min{-^)} envelops. 

2. (cf.p]) The non triviality and stability of the subset $o under the map- 
ping A^* is established in terms of the new basic sequences {6n,max{^)} , 
{5?i,mm(A)} and {Hq} The latter is furthermore used for the proof of 
local contractivity. This proof is simpler in comparison with our pre- 
vious analogs, due also to the fact that we introduce a new norm on 
the Banach space B. 

3. In the present paper, starting from these particular sequences, {Sn,max{-^)} 
and {(5n,mm(A)}, we define a new $o-iteration and explore the behav- 
ior of the Green's functions (essentially the 5 functions), at sufficiently 
large n and reasonable order of this new iteration. 

These last numerical results are convincing, the convergence is rapidly 
established for different values of A and the sign and splitting proper- 
ties {stability + contractivity) give coherent results with respect to our 
theoretical conclusions. 
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2 Reminders 



In [1] we presented in detail the definitions and results introduced in the pre- 
vious papers together with the new ones. Let us present only the necessary 
among them, for the best understanding of our numerical study. 

1 The <1>Q equations of motion the subsets $o C and the new 
mapping Ai* 

Definition 1.1 (The $q equations of motion) 

V A e M+ 

H^{A) = -AH^{A) +1 (1.2.1) 

and for all n > 3, 

//"+i(A) = A"+i(A) + 5"+i(A) + C"+HA) (1.2.2) 

with: 

A"+i(A) = -AF"+=^(A); (1.2.3) 
S"+i(A)= -3A y -^W^+^(A)W^+\Ay, (1.2.4) 

C-^\A) = -6A E J- n (1.2.5) 

Here the notation Wn{J), means the set of different partitions (ji;j2) of n 
such that ji is an odd integer and ji + j2 = n. Respectively Wn{I) is the set 
of triplets of odd numbers-different ordered partitions (ii;i2;^3) of n with 
ii>i2> h- 

The symmetry- integer asym{I) is defined by: 

r 3! if ii=i2= is 

OSymil) = S 1 n / «2 / «3 / n (1.2.6) 

I 2 otherwise 



2 The vector space B 
Definition 1.2 

14^e introduce the vector space B of the sequences H = {H^'^^}n=2k+i;k&n 
by the following: The functions H^^^{A) belong to the space C°°(M^) of con- 
tinuously differentiable numerical functions of the variable A G ( which 
physically represents the coupling constant). 

Moreover, there exists a universal (independent of n and of A) positive 
constant Kq, such that the following uniform bounds are verified: 
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y n = 2k + l,k EN 
\H''+^{A)\ < n !(Ko)" VA e M+ (1.2.7) 

We suppose that the system of equations under consideration, concerns 
always (following our introduction and the previous definition) the sequences 
of Euclidean connected and amputated with respect to the free propagators 
Green's functions (the Schwinger functions), and that these sequences de- 
noted by H = {H"'~^^}n=2k+i ,fceN belong to the above space Be- 

3 The splitting sequences and the subsets $o C $ C B 

Definition 1.3 

We first introduce the class V of sequences 

S = {Sni^)}n=2k+l;ken ^ B, 

such that they verify the bounds (2.1.1) in the following simpler form: 

|5n(A)| < Ko, (1.2.8) 

Vn = 2A; + l; A; e N 

Definition 1.4 splitting and signs-^ C B 

We shall say that a sequence H e B belongs to the subset $ C B if there 
exists an increased associated sequence of positive and bounded functions 

on M+, 

S = {<^n(A)}„=2fc+l;A:eN G V, 

such that the following "splitting" (or factorization) and sign properties are 
verified: 

VAGK+ 



$.1 



$.2 



$.3 



H^{A) = 1 + A(5i(A) with: |im Si{A) = (1.2.9) 



H\A) = -5s{A)[H\A)f, with: dsiA) < 6A, lim = 6 



A-*-0 A 

(1.2.10) 



Vn > 5 H^+\A) = ^-4^")^—.. y^ith : lim = 3n(n - 1) 

^ ^ 3An(n - 1) ' A-^o A ^ ^ 

(1.2.11) 
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$.4 

V n = 2A; + 1 with A; G N*, 3 positive continuous functions of A, 

(A), (A) 

(uniform bounds independent on iJ), ^oo (uniform limit and 

bound at infinity), such that: 

<5n,max(A) > 5n,min{^) (1.2.12) 

'^n,mm(A) < (5„(A) < (A) (1.2.13) 

hm 6n,max < Caa:(A) < (1.2.14) 

Definition 1.5 

S3,maxW = 6A; <^3,m^n(A) = ^ ^ q^^^^ 5^2) (1-2-15) 

and V n > 5 

A <'A^ 3A n(n - 1) n o ifiA 

<5n,max(A) = ■ rv-^ (1.2.16) 

1 + 3A n[n — l)do 



here we put do = 0.001 

n (n — 1) 



3A n(n - 1) n 17^ 

'^"'-^"^^^ = l + 3An(n-l) ^^-^"^^^ 



Definition 1.6 

By using {Sn,max{^)}n=2k+i;km ^"^^ {Sn,min{^)}n=2k+i;ken introduced 
before we define the following sequences: 

{-f^mtx}n=2fc+l;fc6N ^ i3 and {iI^|„^}„=2fe+l;fceN € B 

^^L.(A) = (1 + 6A2)2; i?^,„ = l (1.2.18) 

<ax(A) = -QHhLx?; i?l„(A) = -<53,^,„(A) (1.2.19) 
and recurrently for every n > 5 : 

(A)c;;,+i 

r7-n+l/A\ _ '^n,mm(A)C^"|„ O 01 A 
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Definition 1.7 The subset $o 

Taking into account the sequences of the previous definition we introduce 
the following subset C ^: 

^0 = {HG^:\H;^t^\ < < liJ-tll, yn = 2k + l, k€ n} 

(1.2.22) 

Using the previous sequences we defined in jlj the 'fundamental se- 
quence". 

Definition 1.8 

/f2(A) = l-A<,„ (1.2.23) 

H^{A) = -63,omH'of wtth 53,o(A) = ^AiW^ (1-2-24) 

1 + 9A. I 



\ min I 



and for every n > 5 



r»+l..N_ '^n.o(An"+HA) . 

3An(n - 1) ' 



^n+i(A) = (1.2.25) 



with: 



C^+\A) = -6A Y: . ,. m n ^0 ^'(A); (1-2.26) 

3An(n — 1) , ^ 

and 

7~) / rj \ I mm I I mm I n no\ 

^n[J^{min)) — , rrn+li [L.Z./S) 

\J^max\ 

and proved, 

Proposition 1.1 (the non triviality) 

The set of sequences $ given by the definition \1.4\ is a nontrivial subset 
of the space B. 

Then, we introduced 
Definition 1.9 The new mapping A4* 

M* : 

H^'{A) = 1 + A6[{A) with 5[{A) = -H^{A) (1.2.29) 
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H^'{A) = -5',iA)[H''f wzth 6',{A) = 6A[1 + 6AH\3/2 - g,^^, 

(1.2.30) 

and for every n > 5 

l^n+ll _ 

^n(^) = rjp^: (1.2.33) 



and: 



with 



and proved: 

Theorem 1.1 The stability of the subset $o 

If H e then M*{H) C $ under the condition: 

< A < 0.05 



Furthermore (always in ^), we constructed a Banach space Be <Z B hy 
introducing the following norm M: 

Definition 1.10 

M : B^M+ 

H ^ \\H\\ with: \\H\\ = {svLp- -} ^ ' 

K;n Mn 

Here V A e 

Ml (A) = Hi,, = (1 + 6A2)2; M3(A) = <53max(Mi)3; (1.2.35) 

and for every n > 5 

M„(A)= n{n - l)5nmaxMn-2{Mif (1.2.36) 



The function A/" above defines a finite norm inside a non empty subspace Be 
of This subspace Be obviously contains and is a Banach space with 
respect to the M- topology. 

Using the above definition of the norm we introduced a ball-neighborhood 



Sp{Hq) of the fundamental sequence and by theorem 1.2 , we have established 



the local contractivity of Ai* inside it by the following: 
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Theorem 1.2 the local contractivity of the mapping M* in Sp{Ho) c 

$0 

There exists a finite positive constant A*{k 0.45) such that the mapping 
M* is locally contractive in Sp{Ho) consequently there exiists a unique non 
trivial solution of the $o equations of motion in the neihbourhood Sp{Ho) C 
$0 of the fundamental sequence Hq . 



7 



2 The numerical study 



1 The different aspects of the analysis 

We have studied three different aspects of our study, consequently we have 
obtained three sets of figures, that we describe in detail in what follows. 

The general conclusion of this numerical experience appears clearly the 
same in all three sets. 

We notice that the first three orders of the iteration of 6n,max and 5n,max 
yield different curves which come closer to each other till the fourth order 
iteration. Beyond, i. e. for 5*'* and 6*^* order, we observe a perfect coincidence 

of Sn,max and Sn,min' 

So, when the value of A lies in [0.001, 0.1], the neighborhood where lies 
the fixed point of the contractive mapping is manifestly around the {Hq} 
sequence, {almost first order iteration of {6n,min} sequence). This fact is 
enhanced by the following observation: 

For a given value of A we remark that the sequence {5n,max} decreases 
during the iteration procedure (resp the sequence {6n,min} increases). The 
two sets are almost the same up to the fourth iteration. We notice that the 
decreasing rate of {Sn,rnax} is more important than the increasing rate of 
{<^n,mm}) and this again underlines the fact that the {Hq} neighborhood is 
the best for the local contractivity. 

This result is more satisfactory (from the point of view of the bound of A) 
in comparison with the theoretical proof of the validity of the contractivity 
criteriun at A < 0.05. 

1 First set of figures 

The first set of figures displays the convergence proof of the mapping for 
different values of A, using as stating points both Sn,max and dn,min- 

This 

set represents the results of twenty iterations of the mapping for different 

values of A 

(i.e A G {0.001, 0.01, 0.03, 0.05, 0.075, 0.1}) at n fixed. We have chosen 
ten different values of n: 

n = 7, 9, 11, 13, 15, 17, 19, 21, 23, 25. 

The stability of the values is already attained at the tenth iteration for 
all values of n. 

2 Second set of figures 

The second set displays the convergence of the mapping up to the sixth 
iteration, for different values of A. 
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3 Third set 

The third set of figures displays the summary of the previous configuration, 
for the sixth iteration. 

This figure represents the results of the mapping of 5„ functions, for 

n = 7, 9, 11, 13, 15, 17, 19, 21, 23, 25. 

for different values of A 

(i.e A G {0.001, 0.01, 0.03, 0.05, 0.075, 0.1}) and for all 20 iterations. 

The figure illustrates clearly the convergence of the iteration to the fixed 
point. We remark also that the convergence is more rapid for sufficiently 
small values of A (and even for bigger than the critical point 0.45, due to 
the small values of n). 

Our experience shows that the stability is impossible, for example n = 
1000, when A becomes bigger than 0.05. 

This figure represents the results of the mapping of (5„ functions as sur- 
faces of n and u for fixed A (at the six different values). 

We remark in this figure that: 

• For small values of A (A < 0.001), the "decrease" properties of SnS 
are not apparent. 

• For the intermediate values ("good values") of A, the surfaces show 
the expected concavity as far as u (iteration number) increases. 

• For large values of A (bigger than the critical value ~ 0.045), we 
observe a rapid increasing surface (because we are far from the stability 
and contractivity criteria). 
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3 The figures 



lambda = 0.001 



j.i '•''*..,„..., .... 



Figure 1: 1'''^* set: Convergence up to = 6 with A = 0.001 starting from 

^n,max &nd Sn^min 



lambda - 0.005 



Figure 2: 1^^* set: Convergence up to = 6 with A = 0.0005 starting from 

^n,max and 5ri,min 
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lambda = 01 



sstSIIS 



Figure 3: l^'^* set: Convergence up to = 6 with A = 0.01 starting from 



^n.max S'lld 5^ min 



lambda - 0.015 



Figure 4: 1^'^ set: Convergence up to = 6 with A = 0.015 starting from 



~'n,max 



and 5n 
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lambfia = 0.002 



Figure 5: V set: Convergence up to = 6 with A = 0.02 starting from 



^n.max S'lld 5^ min 



ambda - 0025 



Figure 6: 1^'^ set: Convergence up to = 6 with A = 0.025 starting from 

•>n,inax 



^n.ni.ax and (^^^ mm 
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lambda ^0.003 



°r] 20 40 ■ BO ■ ■ af) 



Figure 7: 1'''^* set: Convergence up to = 6 with A = 0.03 starting from 



^n.max <llld min 



] . tHtmtiihnnttmmtmtmntttmtnm 

f_ °a 20 40 BO ao 



Figure 8: 1'"'^* set: Convergence up to = 6 with A = 0.05 starting from 

^n,max and 5n^min 
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lambda ^0,07 



1 . ; 8 . m » i i 1 1 1 It tut n u u u tmmuunt u tt 
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Figure 9: V^^ set: Convergence up to = 6 with A = 0.07 starting from 



^n.max S'lld 5^ min 



IC 



■ 1 1 1 i i H > » » » 1 1 » > « 1 1 1 » t i i t i i i i t H i i » 1 1 » j 



Figure 10: l^'^ set: Convergence up to = 6 with A = 0.09 starting from 

^n,max and 5n^min 
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lambda - 0.1 



1] °' 



.tUniUtun It I tit II >i u ii itii ii t tit it tt tt » 



Figure 11: 1''''* set: Convergence up to = 6 with A = 0.1 starting from 

^n..m.n.T, and 6ri..nn.in, 



12 



; I ; I « 1 8 N j 8 i I i y 1 1 i 1 1 j I u t.» 1 1 1 1 1.8 ».»i.t « »,» i . i {,i 



Figure 12: l*"^* set: Divergence up to = 6 with A = 0.15 starting from 

^n,max and 5n^min 
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Up to the Sixth delta mapping using Max and Min Lambda-0.001 




Figure 13: 2 set: Convergence up to = 6 with A = 0.001 starting from 



Up to the SiKtli delta mapping using l^aK and Ivlin Lambda-0.01 
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Up to the Sixth delta mapping using Max and Min Lambda=0.02 




Up to the Sixth delta mapping using Max and Min Lambda=0.03 




18 



Up to the Sixth delta mapping using Max and Min Lambda=0.1 




Up to the Sixth delta mapping using Max and Min Lambda=0.15 
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